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ABSTRACT 

We present an analysis of the 101 ks, 2007 Suzaku spectrum of the LINER galaxy 
NGC 1052. The 0.3 — lOkeV continuum is well-modeled by a power-law continuum 
modified by Galactic and intrinsic absorption, and exhibits a soft, thermal emission 
component below 1 keV. Both a narrow core and a broader component of Fe-Ka 
emission are robustly detected at 6.4 keV. While the narrow line is consistent with 
an origin in material distant from the black hole, the broad line is best fit empirically 
by a model that describes fluorescent emission from the inner accretion disk around 
a rapidly rotating black hole. We find no direct evidence for Comptonized reflection 
of the hard X-ray source by the disk above lOkeV, however, which casts doubt on 
the hypothesis that the broad iron line is produced in a standard accretion disk. We 
explore other possible scenarios for producing this spectral feature and conclude that 
the high equivalent width and full width half maximum velocity of the broad iron line 
(v > 0.37c) necessitate an origin within d ~ 8r g of the hard X-ray source. Based on 
the confirmed presence of a strong radio jet in this source, the broad iron line may be 
produced in dense plasma at the base of the jet, implying that emission mechanisms 
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in the centralmost portions of active galactic nuclei are more complex than previously 
thought. 

Subject headings: galaxies:individual - accretion, accretion disks - black holes - 
galaxies rnuclei - X-rays: spectra 


1. Introduction 

Since their first identification as an interesting class at optical wavelengths (Heckman 1980), 
the nature of Low Ionization Nuclear Emission Regions (LINERs) has remained elusive. Many 
such regions are associated with active galactic nuclei (AGN), though their X-ray luminosities are 
highly sub-Eddington and way below those of their Seyfert galaxy counterparts. The reason for 
this modest nuclear energy output remains unknown: possibilities include heavy absorption of the 
nuclear emission, comparably inefficient radiation of the accretion flow, or some combination of 
both. 


The notion of an advection-dominated accretion flow (ADAF) may apply in LINERs, wherein 
the ions lose thermal contact with the electrons and only a small fraction of the dissipated energy 
is radiated (Rees et al. 1982; Narayan & Yi 1994). Much controversy remains as to whether this 
picture is correct (e.g., Guainazzi et al. (2000)), and the debate has become more pointed since the 
report by Ho et al. (1997) that some 20 — 30% of all galaxies are members of the LINER class. 

X-ray spectroscopy can help in determining whether the accretion disk is in an ADAF state 
or not. A LINER galaxy harboring an AGN produces a characteristic hard X-ray continuum. If the 
accretion flow close to the black hole conforms to the conventional thin disk archetype (Shakura & 
Syunyaev 1973), the cool gas orbiting in the inner disk should intercept and reprocess a significant 
fraction of the continuum, also producing the Fe-Ka line (among many fluorescent emission lines 
from other species as well). This spectral feature should have a broadened, skewed profile created 
by the relativistic properties of the spacetime close to the black hole, and a significant Comptonized 
reflection component should be seen in the continuum above lOkeV as well. If, however, the 
accretion flow in this region takes the form of an ADAF, the ion temperature is so high that the 
iron ions should be totally stripped of electrons. Furthermore, the trapped energy will puff up the 
disk, reducing its density and optical depth. Both of these conditions render the reflected emission, 
including the emission of Fe-Ka, virtually non-existent in the inner accretion disk. In this case, 
only contributions to the Fe-Ka line from farther out (e.g., from the outer disk, broad line region or 
putative molecular torus) will manifest, and the line will have a much narrower profile as a result. 

The LINER galaxy NGC 1052 is an interesting case in which to study the dynamics of the 
accretion flow. The AGN is housed by a nearby elliptical host galaxy with a redshift of z = 0.0049 



-3- 


(Knapp et al. 1978), implying a distance of 20.7 Mpc using WMAP cosmology. The source has a 
well-studied LINER optical spectrum (e.g., Gabel et al. (2000)), and is also reasonably bright in 
the X-ray band for a LINER, with an average 2 — lOkeV flux of F 2 _io ~5x 10~ 12 erg cm~ 2 s~ ! . 
The source is known to be heavily absorbed based on the presence of tLO megamasers in its core 
(Braatz et al. 1994), which lie along the direction of a radio jet also present in the galaxy (Claussen 
et al. 1998). A prominent, modestly broad Fe-Ka line has been observed in the source by ASCA 
(Weaver et al. 1999), XMM-Newton, Chandra (Kadler et al. 2004b) and BeppoSAX. This spectral 
feature, coupled with detailed sub-parsec VLBI observations of the jet (Kadler et al. 2004a,b), may 
offer a unique opportunity to study the connection between the accretion flow and jet production 
in AGN. 

A link between accretion and jet activity in AGN has been noted in the radio galaxies 3C 120 
(Kataoka et al. 2007; Marscher 2006; Marscher et al. 2002) and 3C 111 (Marscher 2006). While 
Marscher et al. have noted dips in the X-ray light curves of these sources correlated with ejections 
of superluminal knots in the radio jet, suggesting that a portion of the inner disk is disturbed during 
an ejection event, these results fall short of mapping out the physical structure of the disk at the 
time. However, Kataoka et al. (2007) observed a broad Fe-Ka line in 3C 120 with Suzaku in both 
the high and low flux states of the source. This discovery demonstrates the viability of studying 
the accretion flow and the jet activity of radio-loud AGN simultaneously, and will hopefully pave 
the way for new insights into the fundamental link between these two processes. 

Here we discuss the analysis of our 2007 Suzaku observation of NGC 1052, which offers the 
highest quality X-ray spectrum of this galaxy to date and gives us an unprecedented view of the 
0.5 — 35 keV energy range. NGC 1052 is one of a small handful of AGN sufficiently bright in both 
the X-ray and radio bands for tests of a jet/disk connection. Additionally, it is the only source for 
which we will have a full package of multi-epoch X-ray spectral data, X-ray and radio light curves, 
and VLBA imaging within a common time window and with an analysis by a single collaborative 
group. A forthcoming paper will provide full details of this ongoing campaign, placing our Suzaku 
results in the context of prior X-ray and ongoing radio analysis of this source (Kadler et al. , in 
prep.). 

The X-ray spectrum of NGC 1052 is quite flat and appears to be dominated at lower energies 
by extended thermal and/or photoionized emission (Weaver et al. 1999; Kadler et al. 2004a). The 
source of this emission is likely the interaction of the jet with the surrounding ISM (Kadler et al. 
2004a). Radio continuum and RXTE monitoring both show that the source varies strongly on 
time scales of weeks to months. The 2 — lOkeV flux, in particular, varies between ~ 4-9 x 
10 12 erg cm~ 2 s 1 (Kadler et al. , in prep.), and though there are large uncertainties in the RXTE 
data, the spectrum appears to switch from soft to hard states, perhaps analogous to the high/soft 
vs. low/hard states observed in Galactic microquasars and black hole binaries (Fender et al. 2005). 
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Our Suzaku observation of NGC 1052 and the data reduction will be detailed in Section 2. 
Our approach to the spectral fitting of the time-averaged data and the results of this analysis for 
both the XIS and the HXD/PIN instruments are presented in §3. We examine the Fe-Ka region at 
length in §3 as well and discuss variability over the course of the observation in §4. Our results 
and conclusions comprise §5. 


2. Observations and Data Reduction 

We observed NGC 1052 for a total of ~ 101 ks from 16-18 July 2007 with Suzaku. Reduc- 
tion of the Suzaku data followed the procedures outlined in §4.7 and §4.8 of the Suzaku ABC 
Guide, available online 1 . Calibration files used were the latest version as of the time of this writ- 
ing (September 2008). This calibration was incorporated into the X-ray Imaging Spectrometer 
(XIS) data reduction using the xispi task on the unfiltered data. Once cleaned, the event files 
for the XIS0, XIS1 (XIS2 became defunct as of Nov. 2006) and XIS3 detectors were loaded into 
XSELECT for reprocessing. Images were extracted that enabled identification of source (a circle 
~ 250" in radius) and background regions. After spatial filtering was applied, spectra and light 
curves were extracted for each of the remaining three detectors. Response matrices and ancillary 
response files were generated using the xisresp task, and the spectra, backgrounds and responses 
for the front-illuminated chips (XIS0 and XIS3) were co-added to increase the signal-to-noise us- 
ing addascaspec. The spectra, backgrounds and responses were then rebinned by a factor of eight 
(512 channels) using the rbnpha and rbnrmf tasks to speed spectral fitting. Finally, the spectra, 
backgrounds and responses were grouped together using the grppha task. Only spectral chan- 
nels with with a minimum of 20 cts/bin were used for fitting in order to ensure the validity of the 
chi-squared statistic. This corresponded to an energy range of 0.7 — 9keV for the XIS0+3 data, 
0.5—7 keV for the XIS1 data and 12 — 35 keV for the HXD/PIN data. 

Suzaku did observe NGC 1052 with the silicon diode PIN instrument of the Hard X-ray De- 
tector (HXD), but the GSO crystal scintillator instrument did not have enough counts to provide 
useful data. The PIN data were reduced using the “tuned” calibration files for the non-X-ray 
background (NXB) and the response files that had been generated for the period surrounding the 
observation. These files were downloaded from the Suzaku Guest Observer Facility (GOF) for 
use. We first needed to ensure that only the portion of the NXB coincident with our observation 
was used in processing, so we merged the good-time interval (GTI) of the NXB with that of our 
screened event file to yield a common GTI using mgtime. With XSELECT, we then filtered the 
data and background using this common GTI and extracted our spectra for the source and NXB. 


1 http://heasarc/docs/suzaku/analysis/abc/nodel.html 
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We corrected for dead time in the observation using the hxddtcor task, and increased the expo- 
sure time of the NXB data by a factor of ten to compensate for the inflated event rate, which had 
been similarly multiplied in an effort to suppress Poisson noise. After reduction and filtering, the 
combined XIS+PIN spectrum had 22172 cts from 0.3 — 60keV. 

The cosmic X-ray background (CXB) has not been taken into account in PIN observations, 
so to model it appropriately, we simulated the CXB spectrum in XSPEC (Arnaud 1996) using a 
power-law model of T = 1.29 with a high energy cutoff of foldE= 40keV (initial parameters are 
given in the ABC Guide). We then folded this model through the response file for the PIN flat field 
for an exposure time of 10 6 s to generate the simulated CXB data. This spectrum was then read 
back into XSPEC with the PIN response for our observation and the original model was refitted over 
12 — 35 keV to yield the actual parameters of the CXB during this time. In our case, the simulated 
CXB spectrum contributed a count rate of ~ 2.28 x 10 2 cts s 1 to the X-ray background from 
12 — 35 keV. In modeling the PIN data with XSPEC, we applied the best-fit CXB parameters as a 
constant with each fit. The CXB followed a power-law form with a high-energy cutoff: T = 1.30, 
Apo = 8.13 x 10~ 3 4 phcm~ 2 s _1 , foldE = 40.53 keV. 


3. Spectral Analysis 

We began our analysis of the 2007 Suzaku spectrum of NGC 1052 by examining the data from 
the operational XIS detectors. We use the 0.5 — 9.0keV data to extract information about the un- 
derlying continuum, any soft excess emission, the Fe-K line complex and any intrinsic absorption 
in the system. NGC 1052 is thought to exhibit evidence for an AGN disk (the broad Fe-Ka line) 
and jets (radio observations), as well as substantial, extended soft emission. As such, we expected 
the spectrum to be a conglomerate of all of these features since Suzaku lacks the spatial resolution 
necessary to physically separate them. 

We also analyzed the HXD/PIN spectrum from 12 — 35 keV in order to leam more about the 
nature of the continuum and to constrain the amount of reflection seen in the source. Because a 
broad iron line has been observed in NGC 1052 in the past (Ros et al. 2007), we expected that 
residual emission above a power-law should remain at high energies due to the presence of the 
so-called “Compton Hump.” This feature has been noted in many AGN with broad iron lines (for a 
review, see (Reynolds & Nowak 2003)), representing Compton down-scattering of the hard X-ray 
photons by the material in the disk. This Compton hump typically peaks at ~ 20 — 30keV, and is 
thought to go hand-in-hand with the presence of a broad iron line, as both spectral features arise 
from the same physical process of reflection. 



3.1. XIS Continuum 


Before the iron line region could be considered, it was necessary to first accurately model the 
underlying continuum of NGC 1052. We began by ignoring the energy range where the iron line 
is thought to be important (3 — 7 keV) and modeled the 0.5 — 3 and 7 — 9 keV energy bands with a 
simple power-law modified by Galactic photoabsorption (Ah = 2.83 x 10 20 cm~ 2 , as per Kalberla 
et al. (2005)). The result was a poor fit: % 2 /clof = 840/361 (2.33), where “dof” indicates the 
number of degrees of freedom of the model. Large residuals remained at both high and low energies 
(see Fig. 1), and a very hard photon index for the power-law component was seen: T = 0.34. To 
account for differences in calibration between the XIS0+3 and the XIS 1 detectors, we multiplied 
each data set by a constant for normalization. This constant was held at 1.00 for XIS0+3, and 
achieved a best-fit value of 0.94 for the XIS1 data. 

Especially noteworthy is an apparent soft excess clearly visible below ~ 1.5 keV (see Fig. 1). 
To mitigate this residual feature we incorporated a thermal mekal component (Mewe et al. 1985). 
The fit improved dramatically: % 2 /dof — 381/358(1.06), with kT ~ 0.6keV and a solar abun- 
dance for the mekal component. This component has been noted in previous X-ray analyses of 
NGC 1052 by Kadler et al. (2004b) with Chandra and Weaver et al. (1999) with ASCA (hereafter 
W99), and the temperature and flux obtained from our Suzaku fit are comparable to both values in 
both previous works. This thermal emission most likely originates in extended emission around the 
nucleus from the interaction of the jets with the ISM (Kadler et al. 2004b). Residuals still remained 
at low energies, however, so we added in a second mekal component (with an abundance equal to 
that of the first) to try to account for these features. Although marginal improvement in the overall 
goodness-of-fit was seen (A% 2 /Adof = — 12/2), the addition of this second thermal component ren- 
dered the parameters of both mekal components unconstrained in an error analysis. We therefore 
removed the second thermal component. Though modest residuals remained at select energies, 
no statistically significant improvement in fit was achieved with the addition of discrete Gaussian 
components, as one might expect in a photoionized plasma, or with allowing the abundance of the 
mekal component(s) to vary. Alternatively, we also tried to account for this excess soft emission 
with a second power-law, with and without further absorption. Both models exhibited significantly 
poorer fits than that of the single mekal component. 

Because NGC 1052 is also classified as a Seyfert 2 galaxy, it is reasonable to assume that 
some intrinsic absorption exists in or around the nuclear region. Such absorption could play a 
significant role in the spectral curvature seen, so it must be properly accounted for in the continuum 
model. We employed a partial-covering absorber (zpcfabs) intrinsic to the source to model this 
component, applying it only to the power-law and other components of nuclear AGN emission and 
not to the extended mekal components. A statistically significant improvement in fit was seen: 
% 2 /dof = 335/356 (0.94), with Ah = 1.08 x 10 23 cm~ 2 and a covering fraction of 84%. With the 
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Fig. 1. — The XIS spectrum of NGC 1052 fit with a power-law modified by Galactic photoabsoiption over 
the 0.5 — 3 and 7 — 9keV bands uninfluenced by the fluorescent Fe-Ka line. The fit is poor: % 2 / dof = 
840/361 (2.33). Note especially the large residuals at both the hard and soft ends. Black crosses represent 
the combined XIS0+3 spectrum while the black line represents the model. Red crosses and line represent 
the XIS1 data and model. The lower panel of the figure shows the data-to-model residuals, with the green 
line representing a perfect fit with a data-to-model ratio of one. 
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inclusion of this intrinsic absorption, the power-law photon index increased to T = 0.80. This was 
a significant increase in photon index from the fit without intrinsic absorption. 

Physically, this continuum model and its parameters are roughly comparable to those which 
were fit for the 1996 ASCA observation of the source (W99), although the photon index of the 
power-law is considerably harder and the intrinsic absorption less dense in this 2007 observation. 
A model with a dual neutral absorber modifying intrinsic plus scattered power-law components as 
in W99 did not improve the fit, statistically; moreover, the parameter values of the second power- 
law and absorber could not be constrained. 

Even after modeling the neutral absorption in the system, some residual spectral curvature still 
remained around 2 — 3 keV, indicating that some further component of absorption likely remained 
unmodeled. We applied the simple ionized absorber model absori (Zdziarski et al. 1995) to the 
nuclear emission and noted a marked visual improvement to the fit from 2 — 3 keV, as well as 
a small but significant improvement in the global goodness-of-fit to % 2 /dof = 325/354(0.92). 
This absorber was reasonably thick and moderately ionized, with a column density of Ah = 
1.37 x 10 22 cm 2 and an ionization parameter of c, = 68 erg cm 1 s 1 . With the addition of this 
component, the power-law photon index increased to T = 1.50. This value lies within reasonable, 
expected physical limits for the X-ray continuum of an AGN. Our base continuum model there- 
fore included a power-law modified by intrinsic absorption from both a cold, patchy absorber and 
a moderately ionized absorber of smaller column density, along with extended thermal emission. 
All model components were affected by Galactic photoabsorption. 


3.2. The Fe-Ka Line 

Having fit the continuum, we then took the energies from 3—7 keV back into consideration 
and attempted to model the prominent residuals remaining, highlighted by a strong emission fea- 
ture centered at 6.4 keV in the rest frame that was assumed to be neutral Fe-Ka (see Fig. 2, top). We 
began by holding the continuum parameters constant, except for normalizations and the power-law 
photon index, and refitting. The strong residual at 6.4 keV remained, appearing to have a broadened 
red wing associated with it that extended down to fN j 4 keV. For this fit, % 2 /dof = 711/631 (1.13). 

To account for the strong emission feature of 6.4 keV, we initially assumed some core con- 
tribution from neutral Fe-Ka outside of the inner accretion disk in distant material. This would 
produce a relatively narrow emission line, being outside the region where relativistic broadening 
becomes important. Inserting such a line (a = 5 eV, intentionally less than the resolution of the 
XIS in order to make the line truly narrow) into the model yielded a vast improvement in fit: 
% 2 /dof = 584/629 (0.93), withF = 6.40 ±0.01 keV, indicative of neutral iron, and EW = 111 eV. 
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Fig. 2. — Top: The residual emission feature remaining after fitting the continuum of the NGC 1052 
Suzaku/XIS data, as per §3.1. Note the prominent peak centered at 6.4 keV in the rest frame and the excess 
emission extending down to ~ 4 keV. % 2 /dof = 711/63 1(1. 13). Middle: Residuals remaining after fitting 
only the narrow Fe-Ka core. Evidence for a broader feature is clear. % 2 /dof = 584/629(0.93). Bottom: 
Residuals remaining after a laor line (Laor 1991) was used to model the broad Fe-Ka emission line as well 
as the narrow core. % 2 /dof = 761 /792 (0.96). 


The broad residuals centered on 6.4 keV still remained, however (see Fig. 2, middle). 

We attempted to model these broad residuals three separate ways: with a broad Gaussian, 
with a diskline component representing emission from the inner disk around a Schwarzschild 
black hole (Fabian et al. 1989), and with a laor component representing emission from the inner 
disk around a maximally-spinning Kerr black hole (Laor 1991). Each line is centered at 6.4keV 
in the rest frame, corresponding to neutral Fe-Ka. The Gaussian component improved the fit to 
% 2 /dof = 571/627 (0.91), and was quite broad: o = 747 eV. Its equivalent width of £Wb r0 ad = 
201 eV was slightly smaller than that of previous observations, e.g., W99, in which £Wb r0 ad ~ 
300 eV. Relaxing the energy constraint of the line resulted in E = 5.41 — 5.88 keV, with a = 
0.78 — 1.30keV and .EWbroad = 377 eV, though significant residuals still remained on the red wing 
of the line. Because this is a purely phenomenological model, we turned to diskline and laor in 
the hope of achieving more physical parameter constraints. 

The diskline model yielded % 2 /dof = 574/625 (0.92), which is not significantly different 
from the Gaussian fit, statistically, though it did account for more of the residuals below 6keV. 
VLBI observations of the radio jet on sub-parsec scales have constrained the inclination angle of 
the accretion disk to be between 57 — 72° relative to our line of sight (Kadler et al. 2004b), so we 
applied these constraints to the model. Our best fit yielded an accretion disk emissivity index of 
a = 1.69 — 3.69 (where the disk emissivity £ r~ a ) and an inner disk radius of r m < 26 r g (where 
r g = GM/c 2 ), but the inclination angle could not be further constrained. We held the outer radius 
constant at 400 r g , since it was unable to be constrained and our emissivity index dictated that this 
radius safely encompasses all of the relevant emission from the disk. 

The laor model provided the best statistical fit, at % 2 /dof = 561/625 (0.90), and also ap- 
peared to be the most effective at reducing the residuals between 4 — 7keV (see Fig. 2, bottom). 
Keeping the same disk inclination constraints mentioned above for the diskline model, we again 
could not further constrain this parameter, but we noted that the emissivity index of the disk was 
fairly centrally concentrated at a = 1.53 — 5.10, while the inner radius of emission is only mildly 
constrained within r m = 8.93 — 45.10 r g . The equivalent width of the broad iron line here was 
E VFbroad 185 eV. Interestingly, when we relaxed the constraint on the inclination angle of the 
disk we obtained a best-fit value of i = 45 ± 5°: this is roughly a ~ 2o deviation from the radio 
constraints. 

The laor fit was a slight improvement over both the Gaussian and diskline models, though 
it should be noted that the uncertainty of the r m parameter effectively renders it impossible to 
meaningfully distinguish between the diskline and laor models based on the Fe-Ka profile 
alone. If the laor model fit could better constrain the inner disk radius to r m < 6r g , it would 
suggest that the black hole at the center of NGC 1052 is rapidly rotating, as we would expect for 
a source that powers relativistic jets. We cannot make this statement unequivocally based on our 
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data, however. The data fit with this laor model are shown in Fig. 3. The final best-fit model 
parameters, incorporating the HXD/PIN data as well, are presented in the following Section within 
Table 1 as Model 1 . 

Absorption has often been suggested as an alternative to broad iron emission when attempting 
to model enhanced spectral curvature in the 4 — 7 keV range (e.g., Kinkhabwala (2003)). To address 
this possibility, we compared our best laor fit above to a model in which we removed the broad 
iron line and replaced it with a second absorbed power-law component. The resulting fit yielded 
% 2 /dof = 582/627 (0.93), as opposed to the % 2 /dof = 561/625 (0.90) achieved with the final laor 
model, and left significant residuals around 6.4 keV and below. Moreover, the parameter values 
for both the partial-covering absorber and the power-law components are not constrained in the fit. 
The models employing relativistic disk emission lines unquestionably provide better statistical and 
visual fits to the data. 


3.3. Addition of the HXD/PIN Data 

Adding the HXD/PIN data on to the combined XIS0+3 and XIS 1 spectra, we noticed that the 
signal-to-noise ratio deteriorated significantly above ~ 35 keV. As such, we restricted our energy 
range for this detector to 12 — 35 keV, applying a constant normalization factor as we did with the 
XIS1 data set. In this case, we used a parameter value of 1.09 in keeping with the value expected 
for the HXD nominal pointing position. 

We considered the combined data refitted with Model 1, including no reflection component 
above 9keV. Fig. 4 shows the data, model and residuals from this fit, which yields % 2 /dof = 
660/685 (0.96). The inclusion of the high-energy data left the 2 — 9keV parameters effectively 
unchanged from their values in the fit without the PIN data. If a significant flux contribution 
from reflection is present, as one might assume if the broad iron line has an origin in a relativistic 
accretion disk (see discussion in §3), we would expect to see some excess in the emission over 
9keV relative to Model 1. No such feature was detected, however. 

Table 1 shows the full 0.5 — 35 keV fit with 90% confidence errors. 


3.4. Alternative Continuum Models 

Though no evidence for Comptonized reflection is apparent to the eye above 9keV, it is 
nonetheless important to try to place upper limits on its possible contribution to the overall spec- 
trum. With this in mind, we utilized two popular, well-regarded models of reflection from an 
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Fig. 3. — The best-fitting model for the XIS data. This model includes a power-law continuum modified 
by intrinsic absorption from a partial-covering neutral material as well as a mildly ionized gas, extended 
thermal emission, and both broad and narrow components of neutral Fe-Ka. The spectrum is also modified 
by Galactic photoabsorption. % 2 /dof = 561/625 (0.90). The red, black and green colors/lines are the same 
as in Fig. 1. 
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Fig. 4. — HXD/PIN data included from 12 — 35 keV and refitted with Model 1, which does not include any 
type of continuum disk reflection. Inclusion of this data does not significantly alter the model parameter 
values from the 0.5 — 9 keV fit. Though a broad line is still required at 6.4 keV, no accompanying evidence 
of reflection from the disk above 9keV is seen. Red and black colors follow the designations of Fig. 1. 
In the top panel, the green crosses and line represent the PIN data and model, respectively. In the bottom 
panel the horizontal green line still represents a data-to-model ratio of one, while the green crosses show the 
residuals between the PIN data and model. % 2 /dof = 660/685 (0.96). 
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accretion disk around a black hole in an attempt to quantify the amount of reflection that might be 
present in the spectrum of NGC 1052. 

The pexrav model of Magdziarz & Zdziarski (1995) is one of the most widely used public 
models for AGN emission which incorporates reflected continuum from a neutral accretion disk. 
While it represents a more physical model than the simple power-law alone, it does not include 
line emission from the disk (e.g., Fe-Ka) and it is difficult, in practice, to constrain the three most 
important parameters simultaneously: photon index, reflection fraction and disk inclination angle. 

Adding a pexrav component to our power-law continuum (and keeping the inclination angle 
of the disk constant at the laor value, the power-law cutoff at 100 keV, and the abundances fixed 
at their solar values), we found % 2 /dof = 660/683 (0.97), only a very marginal change in fit from 
that of Model 1 over the 0.5 — 35 keV energy range. The photon index fit between T = 1 .45 — 1 .56 
while the reflection fraction was constrained to be quite small: R < 0.01. As we expect, based on 
the residuals from fitting Model 1, reflection is not a significant contributor to the overall observed 
spectrum. The emissivity index of the laor line for Fe-Ka had very similar constraints to those 
found in Model 1: a = 1.74 — 4.56. Also, to 90% confidence, the inner radius of emission in the 
disk was constrained to r m < 45 r g , as we found in Model 1. The equivalent width of the broad 
Fe-Ka line in this model also remained consistent with that found in Model 1, at EW = 187 eV. 
The pexrav model will hereafter be known as Model 2. 


The ref lion model of Ross & Fabian (2005) is another widely available code used to describe 
the reflected spectrum of an accretion disk around a black hole. The advantage of this model over 
the pexrav model is that it is self-consistent: fluorescent emission lines from many ionized species 
are included in addition to the continuum, most notably those of the Fe-K complex. The ionization 
parameter, incident power-law photon index and iron abundance are specified by the user. The 
disadvantage is that photon indices below T = 1 are outside the allowed parameter range of the 
model, and there is not a parameter for the reflection fraction, as in pexrav. The user is therefore 
obliged to estimate it using, e.g., 


Rrefl 


A|-ef] 7/efl / C, \ 

A po Fp 0 y 30 ) 


( 1 ) 


Here k' re n and K po denote the normalizations of the reflion and power-law components from 
our best-fit spectral model, c, is the best-fit ionization parameter characterizing the disk reflec- 
tion (where t, — 30 denotes a neutral disk), and F re fl and F po are the total fluxes contained in 
the reflion and power-law components for unity normalization over the full wavelength range 
(0.001- 1000 keV). 


We began by removing the laor line from the Model 1 fit and adding in a reflion compo- 
nent, with the iron abundance held at the solar value and the photon index linked to that of the 
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power-law. As expected, this fit mimicked that of the continuum plus a narrow 6.4 keV Gaus- 
sian, leaving residuals on the low-energy end of the line. Statistically, x 2 /dof = 679/687 (0.99). 
Allowing the iron abundance to fit freely yielded a small improvement in fit, with % 2 /dof — 
667/686(0.97) with Fe/solar < 3.53, but doing so rendered the ionization parameter and nor- 
malization unconstrained. As such we elected to keep the iron abundance constant at the solar 
value. 

If NGC 1052 harbors a rapidly spinning black hole, as our best laor fit seems to indicate, 
we must take into account the relativistic effects on emission originating from the inner accretion 
disk. We incorporated this smearing of spectral features using the laor-based convolution model 
kdblur. Our best fit was % 2 /dof = 660/684(0.96), comparable to our best laor and pexrav 
fits. However, we were unable to constrain the emissivity index of the accretion disk, though we 
did obtain constraints on the inner radius of emission that are consistent with those found in the 
laor model: r m < 45 r g . The estimated reflection fraction is quite low at R k q < 0.006, consistent 
with the results of the pexrav model fit and partly attributable to the moderately high ionization 
parameter (£ = 109) in the ref lion model. The kdblur (ref lion) model will hereafter be known 
as Model 3. For comparison, the data-to-model ratios and model components for the Models 1-3 
fits are plotted in Fig. 5, and their best-fit parameters are listed in Table 1. 

Models 2-3 constrain the amount of reflection from the disk to be R ie fi < 0.01, consistent with 
the lack of a discernible Compton hump above lOkeV in the HXD/PIN data. Moreover, as we 
describe in detail later in §5, if the broad iron line results from fluorescence in the disk caused 
by incident hard X-ray emission, then the equivalent width of the Fe-Ka line should be directly 
related to the normalization of the reflection component in our model as compared to that of the 
power-law. Yet this is not the case in our data. Based on the low normalizations for reflection seen 
in both Model 2 and Model 3, we should expect a broad Fe-Ka equivalent width of EW ~ 80 eV, 
much less than the observed EW ~ 185 eV in both Models 1-2. In other words, our broad iron line 
does not appear to have a correspondingly strong amount of continuum reflection associated with 
it. It is possible that this “missing” reflection is actually present to some greater degree in the data, 
but is perhaps swamped out by the power-law continuum. It is also possible that the reflection is 
simply absent in this source. This scenario presents several intriguing possibilities for physical 
interpretations, which we discuss in §5. 


4. Timing Analysis and Variability 


The spectral results presented thus far have focused solely on the time-averaged data from 
the XIS and HXD/PIN instruments, but because many types of AGN have been observed to vary 
on timescales ranging from hours to years, we address the spectral and temporal variability of 
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Energy (keV) 



Energy (keV) 


Fig. 5. — The relative contributions of each model component in the three relevant fits to the NGC 1052 
XIS and PIN data: laor (top. Model 1), pexrav+laor (middle, Model 2) and kdblurf ref lion) (bottom. 
Model 3). Each model presents its own strengths and weaknesses, as discussed in the text. Models are 
extended to 100 keV to illustrate their projected contributions at higher energies. Red represents the power- 
law (affected by Galactic and intrinsic absorption), green represents the thermal mekal component, Gaussian 
emission lines are in dark blue and reflection features such as the laor line (top), the pexrav component 
(middle) and the kdblurf ref lion) model (bottom) are in light blue. The laor line in Model 2 is in 
magenta. 
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Component 

Parameter Value 

Model 1 

Model 2 

Model 3 

mekal 

kT ( keV) 

Kyt ( ph cm “ s' 1 ) 
FkT (erg cm' 2 s' 1 ) 

HUB 



absori 

r 

Ah (x 10 22 cm' 2 ) 

r(R) 

^( erg cm' 1 s' 1 ) 

,-C.I 

1 50 +UU6 

i ' JU -0.05 
1 ^7+0-27 
LJ '-0.23 

3 x 10 4 
66^ 

150+® 

1 ' JU -0.05 

1 45+0.63 
A +°-0.53 

3x 10 4 

67+26 

07 -33 

zpcfabs 

A h (x 10 22 cm' 2 ) 
%cover 

E&Z&nfnMHi 

m rtw 

BflH 

zpo 

r 

Ar(phcm' 2 s' 1 ) 
Ft (erg cm' 2 s' 1 ) 

mm 

Hi 

1KSS 
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E (keV) 

^narrow ( ph Cm ~ S ) 
^narrow ( erg Cm' 2 S' 1 ) 
A VF n arrow ( eV) 

6 40+ o t)1 
°-+ u -o.oi 

1.18±g;g x 10' 5 
l.OOlJ^xlO' 13 

1 1 1 +18 

6 4Q+0-01 
°'+ u -0.01 

1.19±g;J| x 10 -5 
9.54+JJ 4 x 10' 14 

1 1 1 + 18 

6 40+0-01 
°'+ u -0.01 

1.26±g;}| x 10' 5 

1.23±g;}| x 10' 13 

121+ 16 

1Z1 -16 

laor, kdblur 

®emis 

An (r g ) 

K°) 

Abroad (Ph Cm' 2 S' 1 ) 

Abroad ( erg cm' 2 s' 1 ) 

Abroad (eV) 

lllfiBip!] -iJSSti - 
"S'- " \5s~-, 


0 95+ y ui 
u.yj 0.95 

20 15+ 25 - 36 

zu. ro 20.15 

72+0 

7 —15 

pexrav, reflion 

Fe/solar 

£( erg cm' 1 s” 1 ) 

A ref] 

Arefl (phcm' 2 S' 1 ) 
F refl ( erg cm' 2 s' 1 ) 

— 

1.0 

0.0 

0 oi +0 
u - Wi -o.oi 

7.211? 21 x 10' 8 
4.04lj' 04 x 10' 18 

1.0 

109+ 2 

luy -109 

5.75l° 75 x 10' 3 

L30 -L22 x 10_7 

4.24 x 10' 13 


f/dof 

660/685(0.96) 

660/683(0.97) 

660/684(0.96) 


Table 1: Comparison of our three best-fit models for the 0.5 — 35 keV spectrum of NGC 1052. All compo- 
nents of Models 1-3 are modified by Galactic hydrogen absorption with Ah = 2.83 x 10 2() cnr 2 . “K” de- 
notes the normalization value of a given component. Flux values indicate absorbed flux from 0.5 — 35 keV. 
Redshifts for the model components are held constant at the cosmological value for NGC 1052: z = 0.0049. 
Abundances not listed are held constant at their solar values. The inclination angle of the accretion disk 
used in the laor and kdblur models was constrained to fall within the radio observation uncertainties of 
i = 57 — 12° . All errors listed are at 90% confidence for one interesting parameter. 
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NGC 1052 here. The 2007 Suzaku data are discussed below, as are the 1996 ASCA data of W99 
for comparison. The X-ray results from intervening epochs will be addressed by Kadler et al. (in 
prep.). 


4.1. 2007 Suzaku Data 

No significant (> 3a) changes were seen in the count rates of any of the XIS instruments 
or the PIN instrument over the course of the observation. Fig. 6 depicts the light curves from 
XIS 0, XIS 1 and XIS 3, as well as the combined XIS light curve vs. that of the PIN. While 
small variations in count rate (over ~ 10 4 s timescales) were seen in both data sets, neither varied 
more than a factor of two from its baseline flux throughout the observation. This flux was roughly 
0.06 cts s^ 1 for the PIN and 0.38 cts s -1 for the combined XIS data (0.16 for XIS 0, 0.14 for XIS 3 
and 0.21 for XIS 1). 

As mentioned in the previous Section, the overall 2 — 10 keV flux did not vary significantly be- 
tween Models 1-3, averaging to Fi- w = 5.37 x 10~ 12 ergcm~ 2 s' 1 (absorbed). This was split be- 
tween the continuum and the Fe-K lines, where the continuum flux was F con t = 4.86 x 10~ 12 erg cm~ 2 s 
and that of the combined broad and narrow Fe-Ka lines was Fp e = 5. 10 x 10~ 13 erg cm~ 2 s” 1 . The 
average unabsorbed 2 — 10 keV flux at the distance of NGC 1052, F^-io = 8.96 x 10~ 12 ergcm~ 2 s -1 , 
works out to an intrinsic luminosity of L 2 _io = 4.60 x 10 41 erg s” 1 using WMAP cosmology. 


4.2. Variation Between Epochs 

We considered our best 2007 Suzaku results (Model 1) against those obtained in the 1996 
ASCA observation by W99. Though these authors did not present a light curve for their data, 
they fit a time-averaged 2 — lOkeV flux of F 2 _io ~ 8 x 10~ 12 ergcm~ 2 s” 1 , corrected for a dual 
absorber in the best fit. This differs only slightly from the unabsorbed flux of Fi-\o ~ 8.96 x 
10 12 erg cm~ 2 s 1 we have obtained from our data in the present epoch. The neutral hydrogen 
absorbing column in our 2007 best fit is roughly a factor of three less dense than that of W99: 
whereas we find evidence for one intrinsic neutral absorber of Ah ~ 1.08 x 10 23 cm~ 2 , in their 
best-fitting model W99 employ a two-zone absorption model with Ah(1) = 3 x 10 23 cm~ 2 and 
Mt( 2) = 5 x 10 22 cm 2 . Our model contains a component of ionized absorption that the W99 
model does not, however: the absori model representing this absorbing material has a column 
density of Ah = 1.37 x 10 22 cm~ 2 and an ionization parameter of ^ = 68 erg cm -1 s -1 . Even with 
this additional layer the W99 model contains a higher absorbing column than our 2007 model. 

The heavy absorption in W99 resulted in a power-law photon index split between intrinsic 
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Waveband (keV) 

Model 1 ( erg cm 2 s 1 ) 

Model 2 ( erg cm 2 s 1 ) 

Model 3 ( erg cm 2 s 1 ) 

0.5 -2.0 

3.42 x 10 -13 
3.42 x 10" 12 

3.42 x 10 -13 
3.42 x 10 -12 

3.68 x 10~ 13 
5.31 x 10~ 12 

2.0-10.0 

5.53 x 10 -12 
8.96 x 10" 12 

5.25 x 10 -12 
8.96 x 10 -12 

5.33 x 10 -12 
9.20 x 10 -12 

10.0-60.0 

3.64 x 10 -11 
3.68 x 10~ n 

3.64 x 10~ n 
3.68 xlO- 11 

3.79 x 10- 11 
3.82 x 10 -11 


Table 2: Extrapolated model fluxes for the 2007 Suzaku observation of NGC 1052 in three wavebands. 
Absorbed fluxes are on the top line, unabsorbed fluxes are on the line below. For the softest energies the flux 
from XIS 1 was used, while the XIS 0+3 data was used for the 2—10 keV band. PIN data were used for the 
flux above 10 keV. 


NGC 1052 

Bin time: 5000. s 



Bin time: 5000. s 



Time (s) 

Start Time 14297 11:13:25:184 Stop Time 14299 21:33:25:184 


Fig. 6. — Left: Light curves from the three operational XIS instruments for the 2007 Suzaku observation 
of NGC 1052. Right: Light curve of the combined XIS data (black) as well as that of the PIN (red) over 
the course of the observation. All light curves shown are background-subtracted, including the CXB for the 
PIN data. Note the overall lack of significant variability of the source. 
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and scattered components, where a power-law photon index of T = 1.7 was assumed and 11% 
was scattered, while 70% of the intrinsic flux was absorbed by the larger column and 30% by 
the smaller column. We attempted such a neutral dual-absorber fit (see §3. 1-3.2) and found that 
the parameters of the second absorber could not be constrained, and that a single intrinsic neutral 
absorber and an ionized along with Galactic photoabsorption modifying a simple, albeit fairly hard 
power-law (r ~ 1.5) provided an excellent fit to the continuum. 

W99 also found evidence for a thermal mekal component to explain the soft excess noted be- 
low 2 keV, which was best modeled with kT = 0.53+q 26 keV and a 0. 1 — 2.0 keV flux of F meka i = 
5 . 8^2 3 x 10~ 14 ergcm~ 2 s _1 . This component was also detected by Kadler et al. (2004a) with 
Chandra , and we found a similar thermal component in the Suzaku data. Though the XIS de- 
tectors lose effective area rapidly below 0.3 keV, extrapolating Model 1 down to 0.1 keV yielded 
kT — 0.64 ± 0.04 keV and F me k a i = 9.38+2 fi x 10~ 14 erg cm~ 2 s -1 over the same energy range. 
Our results showed a slight increase in temperature and flux for the mekal component but were 
consistent with the 1996 results within error bars. 

In the favored dual-absorber ASCA model of W99, a broad Fe-Ka line component was not 
required in the fit and hence these authors found no definitive evidence for reflection from an 
optically-thick accretion disk in their data. They did require a narrow Fe-Ka line, however: it was 
centered at 6.37 keV and has an intrinsic width of a = 20 eV, an equivalent width of 270 ± 120 eV 
and a normalization of 1.8 ±0.8 x 10~ 5 phcm~ 2 s -1 . In our Model 1 for the Suzaku data, we 
found that both broad and narrow components were necessary in order to eliminate the residuals 
between 4 — 7 keV. Our narrow line, with centroid energy constrained at 6.4 ± 0.01 keV in the 
rest frame with an intrinsic width held at o = 5 eV, had an equivalent width of 1 1 1 ± 18 eV and 
a normalization of 1.18 ±0.19 x 10~ 5 phcm~ 2 s _1 . The presence of a broad component to the 
Fe-Ka line likely rendered this narrow feature significantly weaker than its 1996 counterpart. 

Our broad Fe-Ka line at F = 6.4 keV was best fit with a laor emission feature originating 
from the accretion disk around a rapidly rotating black hole, though the inner radius of emission 
within the disk was only broadly constrained to r m < 45 r g , as discussed in §3.2. This feature was 
stronger than the narrow Fe-Ka line: its normalization was 1.69 +q ^ x 10 5 phcm~ 2 s -1 and its 
equivalent width was 1 85_ 75 eV. Replacing the laor line with a broad Gaussian, we found that 
o = 1.04 ± 0.26 keV, which translates to a lower-limit FWHM of v ~ 0.37c and implies an origin 
for this feature very close to the black hole. Though the fit of the laor line itself is suggestive of 
a broad, fluorescent Fe-Ka origin in the inner accretion disk, the lack of reflected continuum flux 
above lOkeV argues against this interpretation, as we discuss below and in the following Section. 

W99 considered a reflection-based href 1 model to fit the ASCA spectrum, but were unable 
to obtain meaningful constraints on the reflection parameters without making several assumptions. 
These included fixing the value of the power-law photon index of the central X-ray source to be T = 
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1.7 and also fixing the inclination angle of the disk to i = 60°. The latter was consistent with our 
own angle constraints based on the VLBI observations of NGC 1052, but our Suzciku observation 
indicated a somewhat harder photon index (see Table 1). For an assumed reflection fraction of 
R re Q = 1 and the observed ASCA spectrum, W99 found that only a small fraction of the direct 
emission from the X-ray source was visible, leading them to infer that a thick absorbing column 
in the nucleus (Mi > 3 x 10 23 cm~ 2 ) blocked the X-ray source from view, such that reflected 
emission dominated the spectrum. However, the equivalent width of the Fe-Ka line, EW ~ 40 eV, 
was over an order of magnitude smaller than the same line seen in other heavily obscured AGN. 
To account for this line with a reflection model would require an implausibly small iron abundance 
of Zpe/solar < 0.05. This reflection model was therefore argued against by W99. 

The 2007 Suzaku data concur with this assessment of a lack of reflection: though Model 2 
produced a comparable global goodness-of-fit to Model 1, the reflection fraction was found to be 
quite small at /? re fl < 0.01. Attempting to model any reflection present with Model 3, we were 
unable to constrain the emissivity index of the accretion disk, though we did achieve similar con- 
straints for the equivalent width of the narrow iron line and the photon index of the power-law 
component. Based on this spectral fitting, we conclude that in 2007, NGC 1052 does show evi- 
dence for a broadened iron line; this line is not particularly broad or strong, however, and is not 
accompanied by the expected Compton Hump above lOkeV. This may indicate that the line arises 
outside of the innermost portions of the disk and/or that reflection simply does not play a large role 
in creating the overall X-ray spectrum. 


5. Discussion 

We have performed detailed spectral fits to the XIS and HXD/PIN data from the 101 ks, 2007 
Suzaku observation of the LINER galaxy NGC 1052. Three different spectral models were em- 
ployed: (1) an absorbed power-law continuum with contributions from a soft thermal component 
along with narrow and broad Fe-Ka lines; (2) the same model, but with an added component of 
reflected continuum emission from a neutral accretion disk (pexrav); (3) the same base model, 
again, but with a more advanced, self-consistent reflection model which includes fluorescent emis- 
sion lines from the disk and takes ionization into account (ref lion; see Ross & Fabian (2005) for 
details). In Model 3 we applied relativistic smearing using the kdblur model. Because fluorescent 
lines are included in ref lion, there is no need for an additional broad iron line component as in 
Model 1. We summarize our findings below: 

• The three models are quite similar in their statistical goodness-of-fit, though some of the disk 
reflection parameters in Models 2-3 prove especially difficult to constrain. 
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• The continuum is well described by a component with a power-law spectrum of photon 
index T ~ 1.5, thermal emission with kT ~ 0.64 keV, a neutral intrinsic absorbing column 
of N\\ ~ 10 23 cm -2 and a covering fraction of ~ 84%, and an ionized intrinsic absorber with 
A^h = 1-37 x 10 22 cm~ 2 and £, = 68ergcm _1 s~ ! . These findings are relatively consistent 
with those of the 1996 ASCA observation of W99, though our photon index is harder and our 
absorbing column lower in density in 2007. 

• As in the ASCA data, a narrow Fe-Ka line is present: EW = 111 eV. 

• Strong evidence exists for a broadened Fe-Ka emission component (EW = 185 eV). Though 
the best statistical fit is achieved with a laor line profile, however, the rather broad constraint 
on the inner radius of r m < 45 r g renders it impossible to distinguish between a spinning and 
non-spinning black hole. 

• The PIN instrument did detect hard X-ray emission above lOkeV from NGC 1052. This 
emission is well fit by the power-law component of the 2 — lOkeV continuum and does not 
show any evidence for reflected emission via the Compton hump that is commonly associated 
with other AGN harboring broad iron lines, and expected based on the presence of the broad 
iron line if the two features arise from reflection onto a standard accretion disk. Neither of 
our reflection models found any statistical evidence for a strong reflection component. 

• Though VLBI observations of the jet in NGC 1052 have constrained the orientation angle 
of the disk to be between i = 57 — 72°, we find that removing this constraint in our fitting 
yields a best-fit inclination angle of i = 45 ± 5°, though this does not significantly improve 
the global goodness-of-fit. Nonetheless, this result does mark a 2o deviation from the radio 
results. 

These findings motivate three main questions: (1) what is the physical origin of the broad iron 
line in our 2007 observation; (2) how is it that we detect a broad iron line yet no evidence for disk 
reflection above lOkeV, and (3) why does the laor line inclination, when kept as a free parameter 
in the fit, not agree with the disk inclination constraints established in radio observations of the jet 
in NGC 1052 on sub-parsec scales? 

As it turns out, the answers to these questions may be closely related. 

Broad iron lines are expected to be associated with pronounced Compton humps. Comp- 
tonized reflection originates via irradiation of the optically thick inner accretion disk by the hard 
X-ray source, perhaps in some sort of lamp-post geometry, e.g., the base of a jet (Miniutti et al. 
2007). This same process will also produce fluorescent line emission from the disk, most notably 
the Fe-Ka line, which will experience significant broadening and skewness by virtue of its ori- 
gin from the spacetime proximal to the black hole (Reynolds & Nowak 2003). The absence of a 
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Compton Hump, however, casts doubt on the presence/contribution of reflection from the disk to 
the overall X-ray spectrum. From our pexrav model fit, we can estimate the equivalent width of 
the iron line we should expect if both the Compton hump and the broad line originate from disk 
reflection. Using the ratio of the normalization of the power-law to that of the reflected component 
in our fit, we infer that the iron line equivalent width should be EW ~ 80 eV. This low value is 
highly inconsistent with the observed strength of the broad iron feature. 

Relativistic disk smearing is not the only process by which this spectral line can be broadened, 
however. In what follows we consider a number of other physical scenarios, including: transition 
to an advection-dominated accretion flow (ADAF) in the inner disk, resonant scattering from the 
putative torus region, origin of the line in an outflow located in the broad line region (BLR), and 
Fe-Ka produced by Cerenkov line-like radiation rather than fluorescence. 

One explanation for the lack of an observed reflected continuum could be that the accretion 
disk transitions to a radiatively-inefficient state within some radius, e.g., Narayan & Yi (1994). In 
this ADAF picture, which has often been invoked as a potential explanation for low-luminosity 
sources (especially those with observed jet activity), the disk traps radiation rather than releasing 
it, puffing up to become optically thin and geometrically thick as its temperature and ionization 
state rises. Because of these properties, the gas would not be expected to produce much, if any, 
reflection, either in the form of a Compton hump or discrete emission lines. If a contribution to 
the broad iron line did arise from this region, we would expect it to be highly ionized, yet this is 
not seen in the data: our broad and narrow iron line components both have energies robustly less 
than 6.41 keV, suggesting that they arise from species at or below FeXVlll. So while an ADAF 
may indeed be a viable explanation for the lack of a Compton hump in this source, it is likely not 
a significant contributor to the broadening of the Fe-Ka feature. 

The lack of a Compton hump and the mismatch between the unconstrained laor disk inclina- 
tion and the VLBI value causes us to rule out a relativistic disk origin for the broad iron line. We 
also eliminate a possible origin in an ADAF due to ionization constraints and optical depth issues. 
Investigating other possible sources of iron emission in the central AGN leads us to consider the 
role of resonant scattering in the torus region, which is thought to contribute heavily to the flux of 
the narrow Fe-Ka line. However, both the large equivalent width and the large FWHM velocity 
of the broad line component we observe (both substantially greater than the v 2500 km s 1 and 
EW < 100 eV typically associated with the narrow core, e.g., Yaqoob & Padmanabhan (2004)) 
render this explanation implausible. 

An AGN outflow originating in the BLR is explored as a possible explanation for the similar 
spectrum of NGC 7213 (Bianchi et al. 2008). Here a broad iron line with no corresponding Comp- 
ton hump is also observed, and these authors find that the equivalent width of the Fe-Ka feature 
is consistent with predictions made by Yaqoob et al. (2001) for an origin in the BLR or perhaps 
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a Compton-thin torus. Further, measurements of the velocity broadening of the optical Ha line 
concur with those measured from the broad Fe-Ka line (v = 2500 km s _1 ), suggesting a similar 
origin for the emission. Using the same diagnostic, we estimate that we should expect a broad line 
EW ~ 75 — 100 eV if it originates in an outflow from this region; this underestimates our observed 
equivalent width by nearly a factor of ten using a simple Gaussian line to represent the broad 
feature, and by roughly a factor of two if the better- fitting laor line is used. Also, our estimated 
lower-limit FWHM velocity of the line (v ~ 0.37 c) is more than an order of magnitude greater than 
the polarized optical Ha FWHM measurement of Barth et al. (1999) (v 5000 km s 1 ). This is a 
strong indication that the origin of the broad Fe-Ka line we observe comes from a location closer 
to the black hole than the BLR. Moreover, variability of the iron line on years-long timescales has 
been noted (Kadler et al. 2004c; Ros et al. 2007); this type of variation would not be expected 
from the typical symmetric BLR, as changes in flux due to passing clouds in this region should 
theoretically average out over time. 

Finally, we consider a more exotic possibility: perhaps the broad iron line is produced by 
Cerenkov line-like radiation in the centralmost regions of the AGN. You et al. (2003) put for- 
ward the intriguing notion that in dense regions where the refractive index of the material is large, 
relativistic electrons impinging upon this material can produce Cerenkov radiation in a narrow 
wavelength range very close to the intrinsic atomic wavelength of the material. The combination 
of absorption and emission causes the final emission feature to be redshifted and skewed, making 
it appear very like a relativistic line produced by fluorescence in the inner accretion disk. Further- 
more, because electrons are producing the emission rather than photons, no reflected continuum is 
expected. The cool, dense, iron-rich gas in the disk of NGC 1052 could provide the medium for 
this process, while the relativistic electrons in the corona or base of the jet would make for an ideal 
bombardment population. Further, this physical picture makes for a much more consistent match 
with our estimate of the distance of the broad line region of origin from the hard X-ray source: our 
lower-limit FWHM velocity for the broad Gaussian fit to the line yields v ~ 0.37c, which corre- 
sponds to an upper-limit distance of d ~ 8 r g . This is well inside the BLR, and would instead be 
consistent with a lamp-post model for the corona or the base of the jet in NGC 1052. 

Given our constraint on the distance of the broad iron line emitter from the hard X-ray source, 
iron fluorescence from the base of the jet itself is also a potential candidate for the emission mech- 
anism. However, considering the orientation of the jets observed in NGC 1052 to our line of sight, 
it is unclear whether any fluorescent emission from the jet base could escape the system without 
being absorbed. Correlated radio flux measurements from the base of the jet and X-ray spectra of 
this source must be examined over several epochs, simultaneously whenever possible, in order to 
assess the likelihood of this scenario as well as that of the Cerenkov line-like radiation proposed by 
You et al. (2003). Only through a coordinated study of the jet and the accretion flow in NGC 1052 
can we hope to understand the connection between these two vital physical processes in AGN. 



6. Conclusions 


Our 2007 Suzaku spectrum of the LINER galaxy NGC 1052 is consistent with X-ray spectra 
of this source from previous epochs, with a fairly flat power-law continuum that is heavily ab- 
sorbed. Both intrinsic neutral and ionized columns are detected, along with evidence of Galactic 
photoabsorption. An extended thermal component is also detected, likely due to the interaction of 
the jets with the surrounding ISM. Both narrow and broad iron lines are observed, though interest- 
ingly there is little to no reflection seen in the spectrum above lOkeV: two different models used 
to characterize this feature both require R k r < 0.01. 

We are thus faced with a complex scenario for NGC 1052 in which the broad iron line is 
not consistent with the lack of observed continuum reflection from the disk. While it is possible 
that the reflection has been drowned out in some way by the power-law continuum, we must also 
consider the physical implications if reflection is indeed absent above lOkeV. It could be that we 
are witnessing the transition of the accretion flow to a radiatively-inefficient state at some critical 
radius in the disk, which would eliminate the Compton hump from the spectrum. A broad iron 
line is still produced, however, and the velocity width of the line is consistent with an origin close 
to the black hole. If the inner disk is an ADAF, it is possible that the broad Fe-Ka line is emitted 
from material at the base of the jet, or perhaps from a yet more exotic mechanism such as the 
Cerenkov-like line radiation postulated by You et al. (2003). Further, coordinated observations of 
both the inner jet(s) (radio) and the inner accretion flow (X-ray) in order to solve this puzzle. 
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